AN INVERSE SOURCE PROBLEM IN RADIATIVE 
TRANSFER WITH PARTIAL DATA 



MARK HUBENTHAL 

Abstract. The inverse source problem for the radiative transfer equa- 
tion is considered, with partial data. Here it is shown that under certain 
smoothness conditions on the scattering and absorption coefficients, one 
can recover sources supported in a certain subset of the domain, which 
we call the visible set. Furthermore, it is shown for an open dense set of 
C°° absorption and scattering coefficients that one can recover the part 
of the wave front set of the source that is supported in the microlocally 
visible set, modulo a function in the Sobolev space for k arbitrarily 
large. This is an extension to the full data case, where the complete 
recovery of an arbitrary source has been shown. 



1. Introduction 

We consider a problem relevant to optical molecular imaging (OMI), 
which is a fast-growing research area. In this application, biochemical mark- 
ers can be used to detect the presence of specific molecules or genes, and 
suitably designed markers could potentially identify diseases before pheno- 
typical symptoms even appear. The markers are typically light-emitting 
molecules, such as fluorophores or luminophores. In contrast to Single 
Positron Emission Computed Tomography (SPECT), Positron Emission To- 
mography (PET), or Magnetic Resonance Imaging (MRI), optical markers 
emit low-energy near-infrared photons that are relatively harmless to human 
tissue. Further specifics can be found in the bioengineering literature such 
as P El [12] . 

The inverse problem we consider consists of reconstructing the spatial 
distribution of a radiation source from measurements of photon intensities 
at the boundary of the medium in specific outgoing directions. In many 
applications, the propagation of photons emitted can be modeled as inverse 
source problems of steady-state radiative transfer equations. Once we know 
the optical properties of the underlying medium, the problem of determining 
the source is feasible. It is shown in [22] that under mild assumptions on 
the scattering and absorption parameters of the medium this is possible. 
However, in the partial data case, which will be made more clear shortly, one 
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can only hope to recover information about the singularities of the source. 
In particular, we seek to recover information about the wavefront set of the 
source function. We now describe more precisely the mathematical problem. 

We assume $7 to be a bounded domain in M" with smooth boundary and 
outer unit normal vector iy(x). As in [22], we also assume that is strictly 
convex. The authors of that paper note that this is not an essential assump- 
tion, since for the problem considered one can always push the boundary 
away and make it strictly convex without losing generality. Moreover, we 
assume that the data is given on the boundary of a larger domain ili D fi. 
We remark that this condition is not needed for existence and uniqueness 
of the solution to the forward problem, but it is required for the stability 
result ()2.4p . which is adapted from the proof for the complete data case. 

Consider the radiative transfer equation 

9 ■ V^u(x, 9) + a{x, 9)u{x, 9) - j k{x, 9, 9')u{x, 9') d9' = f{x), 

(1-1) u\Q_sn = 0, 

where the absorption a and the collision kernel k are functions with regu- 
larity specified later, the solution u{x, 9) gives the intensity of photons at x 
moving in the direction 9, and d±SQ is the set of points {x,9) € dCl x S""^ 
such that ±i^{x)-9 > 0. That is, d±Sn is the set of points (x, 9) edQx S"-i 
such that 9 is pointing outward or inward, respectively. The source term / 
will be assumed to depend on x only for our purposes. We also remark that 
equation 11.11 is only applicable at a single frequency, as the parameters a 
and k typically depend highly on frequency. In particular, for high energy 
photons there is a coupling between energy and angle, whereas for photons 
with low energy scattering is not accompanied by an energy change. 
In the case of full data, we have boundary measurements 

(1.2) Xf{x,9) = u\9^sn- 

In [22 j . it is shown that for an open, dense set of absorption and scattering 
coefficients {a, k) G C'^ijl x S"-i) x C'^ijl x S^-^ x S"-i), one can recover 
/ G Lp'iVL) uniquely from boundary measurements X f on all of d-^Sil. To 
set up the case of partial data, first let V C ^-^.SQ, be open and let V V. 
Let XV £ CQ°{d+SCl) be a smooth cutoff function such that xv{x,9) = 1 
for (x, 9) G V and xvix, 9) = for (x, 9) ^ V. The boundary measurements 
are then given by 

(1-3) Xvf{x, 9) = xv{x, 9)u\d^s^i- 

To make notation a bit simpler, \iV = d^SVt (complete data) we will just 
write X, since in this case Xy = X. 

In section [3] we will review the direct problem and some relevant results 
of use in the partial data case. We also establish some results about singular 
integrals that will be needed to prove the main theorem. 
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In section H] we consider the inverse problem with partial data, which 
consists of determining the source term / from measuring Xyf- We also 
compute the normal operator XyXy when the scattering coefficient k = 0. 
Note that when a = k = 0, the operators X and Xy are the standard X-ray 
transforms with full and limited data, respectively. When k = 0, then Xy 
is a weighted X-ray transform. 

Following, in section [5] we establish an injectivity result for / € L'^{Q) 
supported in the visible set assuming analytic a. 

Finally, in section [6] we prove the main theorem that one can recover 
the visible singularities of / with respect to the chosen set V when a € 
C°°{Q X S"~^). Both results are based on the microlocal approach used in 
[8]. Results needed pertaining to singular integral operators are located in 
appendix |Al 

2. Statement of Main Results 

When dealing with the inverse problem, which we will describe in detail in 
section we need to take a larger domain (strictly convex for convenience) 
that compactly contains fi. That is, fix a strictly convex open set ili with 
smooth boundary such that D The strict convexity of l^i ensures 
that the functions t±(x,0) are smooth, where t±{x,9) is the travel time 
from X G J7i to dQi in the direction ±0. In other words 

(2.1) {x + T±{x,e)e,e) Gd±sni. 

We will extend a and k to functions on 0,i with the same regularity. We 
choose and fix this extension as a continuous operator in those spaces. Now 
define Xi : L^(r2i) — t- L^{^-^^SQl) in the same way as for X. From this we 
can look at the restriction of Xi applied to functions / supported in f2 by 
first extending such / as zero on Jli \ Q. Essentially, we are moving the 
observation surface outward a bit and taking measurements on d^li. When 
dealing with the inverse problem, we will usually abuse notation and write 
X instead of Xi, with the understanding that we've already extended the 
domain 17 to 17 1. 

It is proven in (Theorem 2, [22]) that the operator X is injective for 
such an open, dense set of coefficients {cr,k) as in (Theorem 1, [22]) with 
/ € L^(r2), and a stability result is obtained for the normal operator X*X : 
L^{n) L2(0). Here the adjoint X* : L^{d+Sn, dT,) L'^{n x S"-i) 
is defined with respect to the measure dS, which we define shortly. More 
specifically, for an open and dense set of pairs {<T,k) S C7^(J7 x S""-*^) x 
C2(n^ X §^r^C"+i(§^"^)), including a neighborhood of (0,0), we have 
that the conclusions of (Theorem 1, [22]) hold in Jli, that Xi is injective 
on L^(r2), and the stability estimate ||/||l2(q) < C\\XlXif\\fjii^Q_^-^ for a 
constant C > locally uniform in (a, k). 

Before stating the main results of this paper, we need to define the set of 
points such that Xy is injective when restricted to sources supported there. 
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Figure 1 . The visible set is shaded here in the case where V 
consists of a connected open subset of the boundary together 
with ah outgoing directions. 

This set wih clearly depend on V. We also denote by lx,eit) as the line 
starting at x with direction 9. 

Definition 1. We define the visible set C 1^ by 

M = {x en\y9 e S""^ 3{z, 9^) e T/with 9^-9 = 

(2.2) such that l^g± intersects x}. 

It is relatively straightforward to show since V is open, A4 is open as well. 
The proof is left to the reader. Now we can state the first main result, which 
is an injectivity condition adapted from results in [8]. 

Theorem 2.1. Let V € d^SQi be an open set and let Ad be as defined 
above. Let W ^ A4. Then there exists an open and dense set of pairs 

(2.3) {a, k) e c\n X s"-i) X c^n^ X sy^; c"+i(s^'-^)), 

including a neighborhood of (0,0), such that for each {o',k) in that set, the 
direct problem (EiP has a unique solution u € L^(r2i x S"~^) for any f G 
L^(n x S"~^), Xy extends to a bounded operator from L^(r2i x S"^-*^) to 
L\d+SQi,dJ:), and 

(1) the map Xy is injective on L'^{W), 

(2) the following stability estimate holds 

(2-4) WfWl^iQ) < CWXyXvfWH^iQi), 

with a constant C > locally uniform in (a, k) . 

Remark 1 . The proof of uniqueness and stability for the direct problem (jl.ip 
as stated in Theorem l2.1l is essentially the same as the one contained in [22], 
so we will focus on the subtle differences. Furthermore, the proof that Xy = 
XvX extends to a bounded operator from L^(Oi) to L^(9+S'r2i, dS) follows 
from the proof in |22j that X is bounded and the fact that multiplication 
by XV is bounded on L'^{d+S^i,dTi). 

For sources / with more general supports, we hope to be able to recover 
certain covectors in the wavefront set of /. Those covectors (x, ^) G T*J7 
that can be detected will depend on V in the following way: 
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Figure 2. An example where (x,^) is in the microlocany 
visible set M' , given that the source / is the characteristic 
function of the shaded set. 

Definition 2. The microlocally visible set corresponding to partial mea- 
surements on is given by 
(2.5) 

M' := {(x,e) G T*n I 36 € S"~^ such that • ^ = and x*(^,^) / 0}. 
Here Xvi^^^) ^^^^ extension of xv '■ dj^S^li M. to Qi x defined 

by xti^^^) = xvix + T+{x,e)e,e). 

Theorem 2.2. Let I be a positive integer. There exists an open dense set 
Oi of pairs {a,k) G C°°(U x §""!) x C^{n^ x x S^"^) depending on 
I such that given {a,k) G Oi, if {z,S,) G M' , then there exists a function 
V G H'-{Q,) such that 

(2.6) {z,0 i WF{X*yXyf) =^ {z,0 i WF{f-rv). 



3. The Direct Problem 
For notational convenience and to be consistent with convention, we set 

(3.1) To = ^ • V^, Ti = To + (J, T = To + a-K, 

where a denotes the operation of multiplication by a{x, 9), and K is defined 
by 

(3.2) Kf{x,e)= [ k{x,e,e')fix,e')de'. 

If /c = 0, we have that 
(3.3) 

Xf{x,9) = I„f{x,e) := [ E{x + te,9)f{x + te)dt, {x,9)ed+sn, 
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where t±{x,9) is the arrival time defined by {x + t±{x,9)9,6) G d±SO, for 
{x,e) enx §''-^. Here E is defined by 



(3.4) 



E{x,e) =exp(^- a{x + se,9)ds^ 



Note that if cr > depends on x only, then is just the attenuated X-ray 
transform along the line through x in the direction 9. Moreover, in this case 
it is injective and [14] gives an explicit inversion formula. 

In the general case with k ^ 0, it is shown in (Theorem 1, [22]) that the 
direct problem (jl.ip is well-posed even for / depending on x and 9. That 
is, for an open and dense set of pairs 

{a, k) G C\n X S"-i) X C\n^ X S^r^; C"+^(S^"^)), 

including a neighborhood of (0, 0), the direct problem Tu = f with u\Q_sn = 
Ohasaunique solutions G L'^{Q,xS^^^) for any / G L^(i7xS"~-'^) depending 
on both x and 9. Furthermore, the complete data operator X, which is only 
a priori bounded when restricted to sufficiently smooth /, extends to a 
bounded operator 

X : L^{n X s"-i) ^ L^{d+sn,dj:). 

The proof of this relies on using the fact that 

(3.5) [T{^f]{x,9) = j exp(^- a{x + T9,9)dT^ f{x + s9,9)ds, 

as well as Fredholm Theory applied to the resolvent (Id — T^^K)^^. 

In order to solve Tu = /, we observe that Tu = Tiu — Ku = /, and so 
applying Tj~^ to both sides yields 

(3.6) u = T{^{Ku + f). 
This is equivalent to the integral equation 

(3.7) (Id - T^^K)u = T{^f. 

Thus, if Id — T^^K is invertible, we can solve the forward problem uniquely 
for 

(3.8) u = T-^f = (Id - T^^Ky^T^^f. 

To find k such that T~^ exists, we note that {Id — T^^ K^^T^^ = T^^(ld — 
KT^^)~^ and look at the operator 

(3.9) ^(A) = (Id - (Ai^Tf 1)2) 

It is shown in |22) that the operator (^KT^^^^ is compact, and for A = the 
resolvent (13. 9p exists. By the analytic Fredholm theorem (Theorem VI. 14, 
[TBj). we have that ^(A) is a meromorphic family of bounded operators with 
poles contained in a discrete set. It can be shown that 

(3.10) (Id - XKT^^y^ = (Id + XKT^^) A{X). 
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In particular, the r.h.s above is a easily seen to be a right inverse. To show 
that it is a left inverse as well, we can expand ^(A) as a Neumann series for 
||ii'T]^^|| <^ 1 and then use analytic continuation to show that it remains 
true for all A that are not poles of ^(A). These ideas will be useful later 
when proving Theorem 12.11 

4. The Inverse Source Problem with Partial Data 

Let V C d+Sil. be some open subset. Then the boundary measurements 
for the problem (jl.ip with partial data are modelled by 

(4.1) Xvfix, 0) := xv{x, 0)u\9^sn, {x, 6) G d+S^ 

where u{x,6) is a solution of (jl.l|) . and xv ■ S+Sil — )■ [0,1] is a smooth 
function equal to for {x,6) ^ V and xv{x,0) = 1 for {x,6) G V V for 
some open V. We also define the operator J : L^(0) — > L^(S7 x §"^^) by 

jf{x,e) = f{x). 

If A; = 0, we have that 

(4.2) Xvfix, 9) = hyf{x,e) := xv(.x,e)I^f{x,e). 

We will proceed as in [22] by looking at Xy as a perturbation of Ia,v- 
Wishful thinking suggests that XyXy is a relatively compact perturbation 
of I*y/o-,V) the normal operator corresponding to /c = 0. Here Xy is the 
adjoint of Xy with respect to the measure on d^SQ given by 

(4.3) d^ = \e-u{x)\dS-^dSe, 

where as stated earlier, i/{x) is the outward unit normal to the boundary 

an. 

First let us consider the case when A; = and compute I*y- Note that 
Ia,v ■■ L^{^ X S""-^) L^{d+Sn,d^), and hence I*y : L^(d+Sn,dJ:) 
L'^{Q X S"""*^). For now we will restrict ourselves to applying I^-y to functions 
/ that depend on x only. Given h{x,6) G L^(9+5r2, dS) and /(x) G L'^i^), 
one can show that 

{i:yh{x),f{x))L2^nxf^n^^)= [ [ h*{y,9)xt{y,0)E{y,9)7{y)d9dy, 

where g'^{x,9) is the extension of g : ^ M to 17 x S"^-*^ defined by 

g'^{x,9) = g{x + t+{x,9)9,9). We also made use of the diffeomorphism 
(j) : d+Sn xO ^nx S"~i where O = {{t_{x,9),0) \ {x,9) G d+SQ}. This 
map is defined by (l){x,9,t) = [x + t9,9). The Jacobian determinant of 4> 
is \v{x) ■ 9\] see (Lemma 2.1, [5j). Note that (j)'^ : Q. x S"-i ^ d+SVt x O 
is given by (j)~^{x, 9) = {x + r+(x, 9)9, 9, t^{x, 9)). Hence the adjoint in the 
no-scattering case has the equation 

(4.4) Ilyh{x,9)= f E{x,9)h*ix,9)xtix,9)d9. 
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4.1. The Normal Operator I*ylcjy. Similar to the way in which we 
derived the adjoint operator /*y, we may compute the normal operator 



{K,vIa,vf{x),g{x))L2(^Q) 

E{x,e)xt{x,0)(.ia,vf(.x,e))* dd 



g{x) dx 



(4.5) 



E{x, e)xt{x, 9) ( xv{x, 9) I E{x + t9, 9)f{x + t9) dt 

2 



E{x, 



# 



g{x) d9 dx 



\y-x\' 



\y-x\' 



n Jn 



\y -X 



n-l 



-g{x) dy dx. 



In the last line we used the substitution y = x + t9 to convert the integral 
over S"""^ X R in {9,t) to an integral over $7 in y. Thus 
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(4.6) I*yl^yfix) 



E{y, 



y-x - 

\y-x\ ■ 



\y-^\ 



y-x ' 
\y-x\ ' 



f{y) 



|n-l 



dy. 



\y - x\ 

In the case that a is C°°, we would like to know where I* yla,v is elliptic. 



One can show using (Theorem 3.4, pj) that 



(4.7) 

where 
(4.8) 



i:yl^yf{x) = {27r)-^ I , 



,i{^-y)< 



b{x,Ofiy)dydC 



bix,C) = (27r)- 



-i{x-y)-(_ 



E{y,§E^i)E{x, 



y-x ■■ 
y-x\ ' 



# / y—x 



\y -X 



\n-l 



dy. 



We can now see that for a G C°°(0 x S" ^) I*yla,v is a pseudodifferential 
operator of order —1, since its kernel is weakly singular. More specifically, 
Proposition 1 of [22] gives that I*yla,v ■ L^i^) ^ H^i^). See also [23]. 

Now, unfortunately equation ()4.8p isn't particularly useful when trying 
to determine where b{x,^) is elliptic. But recall (j4.5p . which shows that 



(4.9) 



KyIo,vf{x) 



A{x,t,9)f{x + t9)dtd9, 



for a particular function A. By Lemma 4.2 of [8] we have that if ^ G 
C^iQi X M X §"-!) (which occurs if a € C~(Qi x S^^-^)), then /* yl^y is 
a classical ^'DO of order —1 with full symbol 



(4.10) 



m=0 
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where 

■m f 

(4.11) 6„(x, i) = 2^— / dTA{x, 0, ^)(^(") {9 ■ dO. 

ml J^n-l 

To check for elHpticity, we need only look at the principal symbol corre- 
sponding to when m = 0. This is just 

(4.12) boix,0 = 27T [ \Eix,e)\''\xtix,9)\^dSi9). 

Since E is nonvanishing, we immediately have by (j4.12p that b{x, ^) is elliptic 
on the set A4' . 

5. InJECTIVITY of Xy RESTRICTED TO THE VISIBLE SeT 

5.1. Injectivity of /o-y and Since we are only able to access 

some open subset V of d+S^li, we cannot expect for the operator I^y or 
the normal operator /* ylay to be injective. However, from [8] we can 
obtain injectivity for sources / supported in a particular subset of fi. But 
first we must introduce the notion of a regular family of curves. We will use 
the notation Ix.e to denote the line segment through a; G in the direction 
9 € S""^ with endpoints on d^i. We can also assume that lx,e{0) = x and 
^'x ei^) ~ ^- clear that the lines Ix^d depend smoothly on (x, 9) € T$7 

in the sense that the function l{x, 6, t) = Ixfiit) depends smoothly on x, 9 and 
t separately. In fact, we have l{x,S^,t) = x + tS, where t £ (a(a;, ^), 6(x, ^)), 
an interval containing 0, and /(x, ^, a(x, ^)), Z(x, ^, &(x, ^)) G di^i. 

Definition 3. Let T be an open family of smooth (oriented) curves on 0, 
with a fixed parametrization on each one of them, with endpoints on dQ, 
such that for each (x, ^) G TCl\0, there is at most one curve 7^^^ € T through 
X in the direction and the dependence on (x,^) is smooth. We say that 
r is a regular family of curves, if for any (x,^) € T*Q, there exists 7 E F 
through X normal to without conjugate points. 

Remark 2. In our specific case, all curves taken are straight lines and have 
no conjugate points. Moreover, if we let Tj^ be the set of lines in ili which 
intersect A4, then it turns out (as is shown in the proof of Theorem 15. ip 
that r_A4 is a regular family when restricted to A4. This is the motivation 
behind how A4 was defined in the first place. 

Theorem 5.1. Let a be analytic on Cl. If I^yf = for f € 'D'{Q,i) sup- 
ported in W (S: A^*"*, then / = 0. In particular, I^y is injective on L^[W). 

Proof. It is clear that the collection T of lines in ili is an analytic regular 
family of curves (see [8j). Let Tj^ be only those lines which pass through 
A4. We claim that T_m is a regular family of curves when restricted to A4. 
To see this, let x E and 9 E S^~^. Then by definition of Ai there exists 
z G V and an angle 9^ normal to 9 such that the line g± passes through 
X and {z,9-^) E 9+50. By definition we also have that /^ g± E which 



10 MARK HUBENTHAL 

proves the claim. Now, by Theorem 1 of [8] we have that / is analytic on Ai 
with support properly contained in A^'"*. In particular, / = on an open 
subset of each component of A4. Therefore / = 0. □ 

Although the definition of is a bit cryptic and difficult to visualize, 
it is possible to easily visualize an important subset of M when V has a 
certain form, as shown by Lemma 15.21 Here we use the notation chA to 
denote the closed convex hull of a set A C M". 

Lemma 5.2. Suppose that V = 7:~^{W) where W consists of a countable 
collection of disjoint connected open subsets of dQi, and it : d^SO,i — t- d^li 
is the natural projection. Then [jj (chW^)™* C A4 where Wj is a given 
component ofW. 

Proof. Suppose W = |J„ Wa where Wa C dVt, are disjoint connected open 
sets. Let x G (chl^Q,)'"*. Let 9 € 5*^"^ and let 9-^ be any vector perpen- 
dicular to ^. If we consider that 

chWa = {hyperplanes P C | P n W„ = 0}^ 

then l^ qa_ must intersect Wa at some point z. Changing the direction of 9^ if 
necessary and using the strict convexity of ili, we have that (z, 9^) S d+Sili. 
This proves that (ch^Vo)™* C A4 for all a. □ 

5.2. Computing Xy as a perturbation of I^^y for k ^ 0. In order to 
approach the case that fe 7^ 0, we will compute explicitly how Xy differs 
from lay- Note that 

(5.1) Xf = xvR+T-'f = XvR+{ld - T,-'K)-'T^'f, 
where 

R+h = h\Q^sn- 

If / depends on x only (the case we are primarily interested in), then 

(5.2) Xyf = xvR+T-\jf = Xy^+(Id - T^-^K)-%-^Jf. 
Now consider the identity 

(5.3) (Id - T^^K)-^T^^ = Tf ^(Id - iCTf ^)-\ 
which implies that 

(5.4) Xyf = xvR+T^-Hld - KT^^y^Jf. 
Writing Xy = I^y + Ly and noting that 

hyf = XvR+Ti^Jf, 

we have that 

(5.5) Xy = lay + xvR+i-U + (Id - TyK)-^)T^^J. 
and so we have 

(5.6) Ly := XvR+{-U + (Id - T^^K)-^)Tyj. 
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Furthermore, 

(5.7) XyXv = I*yIa,V + -^V) C,v '■= K,v^v + Lyicry + LyLv- 
Lemma 5.3. The operators 

dxI*yLv, dxLylay, d^LyLv 
are compact as operators mapping L'^{Qi) into L'^{Qi). 
Proof. Following the steps of the proof of Lemma 3 in [22j, first note that 

(5.8) (-Id + (Id - Tf ^i^)~i)rf ^ = T^^KT^\ld - KT{^)~^ . 
To prove it, we note that 

T^^KT^^ = (-Id + (Id - T{^Ky^)T^^{ldi - KT{^). 
Thus Ly can be written as 

(5.9) Lv = xvR+T{^KT{\l(i - KT{^)-^J. 

We note that multiplication by xv to obtain Ly from L is bounded 
and hence preserves compactness. First we need to analyze I* yLy = 
I*yXvR+T{^KT{^h, where h = h{x,e). Recah that 

[i:yh]{x)= [ E{x,e)h*ix,9)xtix,9)de. 



Again, as in [22] we notice that xvR+Ti looks like lay, except that now 
the source depends on 9 and x. Thus 

[i:yXvR+T{^9\{x)= I E{x,9)\xv{x,9) f E{x + t9,9)g{x + t9,9)dA d9 
(5.10) 



\y~x\ J y I i'' \y~x\ 



where Fc^cn{x,9) is the even part of F with respect to 9 (i.e. Feven(3^)^) = 
i(F(x, 9) + F{x, -9)) ). To get back to I*yLv, we can let g = KT^^h. 

To proceed, we will now make a slightly weaker assumption on k than 
stated in (j2.3p (see [22]). We will assume that k can be written as the 
infinite sum 

oo 

(5.11) k{x,9,9') = ^Gj(^)Kj(x,e') 

where Qj and kj are functions such that 

oo 

(5-12) ^ ||0ilUi(§"-i)||/^j||L°°(nixS"-i) < oo 

i=i 

In particular, we could take Qj to be the spherical harmonics Yj, and 
then Kj would be the corresponding Fourier coefficients in such a basis. 
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As discussed in [22], uniform convergence of (jS.lip is guaranteed if k € 
L°°($7i X Sg,"^; Cg"*"^(S"~-'^)), which is indeed a weaker assumption. 

Now let Kj be the integral operator with kernel QjKj and Bj = kjT^^, 
where we regard nj as integration in 6' against the kernel nj. Thus, 

(5.13) lK,T^^h]{x,e)=0j{e)iB,h](x). 

,5.14, BM^) ^ I ^) 

By the proof of Lemma 1 in [22], we have that Bj{ld-KT^^y^J : L^i^i) 
L^(r2i) is compact. Now observe that 

oo 

(5.15) = ^[d,i:yXvR+T^'ejJ] [Bjild - KT^Y'J] 

i=i 

By ([51^ and Proposition 1(b) of p], we have that da:I*yXvR+Ti^QjJ ■ 

L^(ili) L^(Qi) is bounded with a norm bounded above by C||(T||^2^Qx§n-i') || ||hi(§"-i) • 

Thus each summand of (j5.15p is a compact operator with norm bounded 

above by C||0j||j|^i ||Kj||2,oo, with C depending on a. By the condition (j5.12p . 

we have that dxl* yLy is compact. 

Now, the proof for dxLyLy is similar. In light of the fact that i?j(Id — 
KT^^)^^ J is compact, it suffices to show that dxLyXv R+Ti^ J is bounded. 
Note that KT^^ commutes with (Id — i^T^^^)^^, and hence 

(5.16) 

L*yXvR+T{^QjJ = {xvR+T^^KT^\ld - KT^^r^jf xvR+T^^&jJ 
= {KT{^JY{xvR+T{\ld - KT^Y'TxvR+T^'J- 

As proven in [22j, by the boundedness of xvR+Ti^, the compactness of 
dxLyXvR+Ti^ J relies on dx{KT^^ J)* , and indeed it is. 

Finally, to show that dLylay is compact, we can proceed similarly to 
the case of dLyLy- Observe that dxLyIa,v = LyXvR+Ti^J, which is 
equivalent to (|5.16p with Qj = 1. □ 

Now we are ready to prove the Theorem 12.11 regarding the injectivity of 
Xy when restricted to sources / supported compactly in the visible set Ai. 

Proof of Theorem \2.1\ Our proof mostly parallels the proof of Theorem 2 in 
[22] . By Lemma [5^3] we have that XyXy is equal to /* yl^y plus a relative 
compact operator Cy- First assume that a and k are C°°. In this case, 
I*ylay is elliptic on and thus there is a parametrix Q of order 1 which 
we view as an operator Q : H^{Cli) — )• L^(Q) (We've restricted the image 
to ri, though Ai would do based on our assumption on the support of /). 
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Thus, for / supported in W A4, we have 

(5.17) Qllyl^yf = f + K^f, 

where Ki is of order —1 near M.. Now apply Q to XyXy to get 

(5.18) QX*yXvf = f + Kif + QCvf=:f + K2f. 

By Lemma [5. 31 we have that QCy is compact. Furthermore, Ki : L?'{M.) — > 
LP'{M.) is compact by Relhch's lemma since it is smoothing near M. This 
reduces the problem of inverting XyXy to a Fredholm equation. By Theo- 
rem lS.H we have that for a real analytic on x (S"~^), I^jy is injective when 
restricted to / supported \n.W M. From this point, the proof follows the 
same as that for Theorem 2 of |22j . and so we conclude. □ 

6. A MiCROLOCAL Result 

Although injectivity is a bit much to ask for in the partial data case, it 
is possible to analyze how singularities are propagated under the normal 
operator I*yl„y. Assuming suitable smoothness conditions on fc, we will 
prove that one can partially recover the wavefront set of /. This is somewhat 
analogous to Proposition 1 in [8], which allows one to partially recover the 
analytic wave front set of / when /c = 0. 

To get more of an intuition for how to proceed, consider the case when 
||rj~"'^i(r|| < 1. Since Xy = xv'-R+Tj~^(Id — KT^^)~^J ^ we may use a Neu- 
mann series expansion to get 

(6.1) Xv = XvR+T{^ [yS^^i^)^ J- 

The term corresponding to j = is exactly Ia,v, which does not account 
for scattering. Subsequent terms in the expansion incorporate scattering of 
higher and higher orders. So if we can show that KT^^ is smoothing, then 
the most singular part of the data will be captured in the ballistic term. 

Before proceeding, we need to establish a bit of notation. We define the 
space T-Liifl x S"~^) as the completion of C°°(0 x S"^^) with respect to the 
norm || • WuiinxS"-^) given by 

OO km,n 

(6.2) ||9(a;,^')llw;(f2xS"-i) •= X] X] l|am^ll/f'(OxS"-i)ll^m'^nllHi(S"-i)' 

m=0 k=l 

where g{x,e) = JZ^^ J^lZi a'ni\x)Yi^,l{e) is the series representation of g 

with respect to the spherical harmonics Ym}i{0) (see appendixEl). Similarly, 
we define the space C;(i7 x S"~^) as the completion of C°°{^1 x S"^^) with 
respect to the norm || • ||c,(rix§"-i) given by 

(6.3) \\9{x,0)\\c^(nxS"-i) = ^ Yl \\^m\\c'{n)\\Y-^l\\HHS"-i)- 

m=0 k=l 



14 MARK HUBENTHAL 

The following lemma establishes the regularizing properties of the opera- 
tor KT^\ 

Lemma 6.1. Let f G Hiini x S"-i) with supp{f) C J7 x S"-i and I > 1, 
and suppose that a G C°°(n x W'-^) and k G C°°(n^ x S^r^ x Sg"^). Then 

Proof. First let us recall from the proof of Lemma 1 in [22] that 
(6.4) 

(^'|2^~y|'^=||)^(^'^'|l=f|) / x-y 



X, \x- 



where s, 9') = exp ^— J^^ a{x + r^', 0') dr^ . The characteristic Si 

y\,9')k{x, 6, |f5|y) satisfies the hypotheses of Proposition lA.il and the result 
follows. □ 

Corollary 6.2. Suppose that a G C°°(n x S"-i) and k G C°°(n^ x S;^r^ x 
S^"^). Then {KT^yjf : ^^(Jl) ^ HjiQi x /or a// j > 0. 

Prom this result, we see that in the case that ||i^T-|^^|| < 1, Xyf is equal 
to Ia,vf plus a remainder consisting of a series of terms with successively 
higher regularity, corresponding to higher order scattering. 

Lemma 6.3. Suppose that a G C~(f^xS"-i) and k G C°°(17xS"-i xS"-i). 
Then KT{^K : H^^i x 8"-^) ^ x S"^i) is compact /or all I > 0. 



Proof. Recall from the proof of Lemma 2 in |22j that 
(6.5) 

r r a (x,y,\x - y\,-^^,e,0') 
[KT,-'Kf]{x,6)= / / ^ -i^^ ^-f{y,e')dydd' 

with some C°° a compactly supported in x and y. The integral in y is a 
weakly singular integral of the form in Proposition lA.H and so by part (a) 
we gain a derivative in x for each fixed 9'. Moreover, the smoothness in 
6 of KT^^K f{x, 9) is dependent only on the smoothness of a. Therefore, 
KT{^K : H\ni X W^) i?'+^(Qi x S"-i). By Rellich's Lemma, the 
inclusion H^~^^{Q.i x S"~^) H\Q.i x S"-"^) is compact, which completes 
the proof. □ 

Again suppose that a and k are C°°. Since Ly = xv R+T^^ KT^^ (Jd — 
KT^^y^J, we have by Lemma O that Ly : H^Q) Ui+iiO.i x S"-i). 
Since XyXy = I*yI(T,v + ^Cy and for smooth a, I*yla,v is a pseudodiffer- 
ential operator of order —1, we would like to show that Cy maps H^{Q) into 
H''~^'^{^i). We have the following proposition: 

Proposition 6.1. Let I be a positive integer. There exists an open dense 
set Oi of pairs (a, k) G C°°(n x §"-i) x C7~(n^ x S^r^ x S"^'^) depending 
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on I such that for all < I' < |, the operator 

(6.6) Cv = KyLv + Lylay + LyLy = Ky^v + LyXy 

maps H^'{n) into H^'+^{Qi). Moreover, we can write Ly = F + R where F 
is a pseudodijjerential operator of order —2, and R : L^(J7) H''{Q). 

Proof. First we write 

^vf = Ky^vf + LyXyf 

= {xvR+T{\jyxvR+T{^KT{\lA - KT{^)-^Jf 

(6.7) 

=: hf + hf. 

Given equation (|5.10p . by Proposition lA.ll we have that (xyiJ+Tj^^ j)* xvR+T^^ 
maps nv+i{9.i x E^-^) into H^'+^iQi) for ah 1' > 0. 

We claim that for an open dense set of (a, fc) G C°° x C°°, {ld-KT^^)-^J 
maps H^'i^i) to nvi^i x §""1) for all < /' < Z + 1. First note that by 
Lemma [6jA(i^rf^)2 : H^' (Qi x S^-^) H^' x S""^) is compact for 
all V > 0. Using the analytic Fredholm theorem on the resolvent A{\) with 
(j3.9p and (j3.10p . we conclude that (Id — Ai^T^"^)^^ exists and is bounded on 
H^'i^i X S"-i) for ah A in some complex neighborhood of [0, 1] except for 
possibly a discrete set, which depends on V . Taking the complement of the 
union of all such discrete sets for </'</ + !, we obtain that (Id— i^T^^^)^^ 
is bounded on each H^' {Vli x S"^^) for all < Z' < / + 1 and for all but a 
discrete set of A. So the set of pairs (cr, k) G C°° x C°° for which (jl.ip has 
a unique solution and (Id - KT{^)-^ : H^' (Qi x S"^i) ^ H^' {^i x S"-i) is 
bounded for < /' < Z + is open and dense. Now, we just apply Lemma [6. II 
to the KT^^ factor in /i, which shows that Ii maps W{VL) into H^'+^{ni) 

for o<r <l. 

To analyze I2, we will use the series expansion of k{x,9,6') previously 
defined in (j5.1ip . Observe that 

^2/ = {f^XvR+Tr'KjT^'ild - KT^Y'A {xvR+T^\ld - KT^^y^j) 



(6.8) 
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Similar to (jS.lOp . we can compute 
(6.9) 

(y^ §E^) 0. («) E (y, ^) 9 (y, ^ 



By Proposition lA.H it is then evident that [xvR+T^^&jJ)* XvR+Ti^ 
maps Hviyi X S"~^) into H^'~^^{Q) for any V > 0. Applying Proposition 
lAHto gives that Bj{ld - KT^^y^J maps H^'+^{^i) to i7''+2(Jli) 

for < I' < I. Therefore, so does its adjoint. Altogether, we have that 
hf G H^'+'^i^i) for / G H^' {n) where < T < L 

For the next part of the proposition, note that for m > 1 

m 

(Id - KT{^) Y,{KT{^yj = J - (i^Tf i)"^+V. 
i=o 

By Lemma l6. II each summand {KT^^y J is a pseudodifferential operator of 
order —j with symbol depending smoothly on the parameter 0. Therefore, 
we may construct a symbol 

oo 

(6.10) p{x,i,e) r^Y.^L{{KT,-yj){x,C,9), 

3=0 

where cjl : L^g(r2) 5'^g(r2 x M") is the full left symbol map. Here we 
are using the notation U^p{Q) to refer to pseudodifferential operators with 
symbols in the class 5*^(0 x M") (see 0). The symbol p corresponds to a 

pseudodifferential operator F of order with smooth parameter 9, and we 
have 

(Id - KT-^) oF = J + Ro, Roe (L-;5° (J]) ; C7°^ (§"-!)) 

(6.11) F -R = {Id- KT^^y^J 

R={ld- KT^y^Ro : L^{n) H\n; C7~(S"-i)). 
Substituting (j6.1ip into the expression (j6.7p for Cy gives that 



■dy. 



jCv = 
(6.12) 



F + R, 



R+T{^jy XvR+T{^KT{^F + {xvR+T{^KT{^fX {xvR+T{^f) 1 + R 



where R involves all the terms with R. It is then immediate that F is a 
pseudodifferential operator of order —2 while R maps L^(il) — )• H\VLi). □ 

For reference, given an operator A S L^q{VL), we define WF(A) as the 
smallest closed cone C C T*0 \ such that o"a|c= € S'~°°(C'^), where a a is 
the symbol of A. We are now ready to prove the main theorem: 
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Proof of Theorem \2.2l Assume first that {a, k) € C°° x C°° is in tlie same 
open, dense set as in Proposition 16.11 Since a is C°°{Q x S"~^), we have 
that I*yla,v is a pseudodifferential operator of order —1. Furthermore, it is 
elliptic on N*l{xo,eo) by ([23]). Let e N*l{xo,0o). Then there exists a 

microlocal parametrix Q G L\q{Q,i) elliptic at (z,^) and Si € L^q{Q,) such 
that 

Q/;^'/,^ = id + 5i, 

and (z,^) ^ WF(5i). We will also restrict the image of Q so that Q : 
H^iili) L'^in). Since WF(5i/) C WF(5i) n WF(/) (e.g. by Lemma 
7.2 of [9]), we have that 5i/ is microlocally smooth near (z,^), i.e. {z,^) ^ 
WF(Si/). 

Now we apply Q to the normal operator XyXy = I* yla,v + to get 
QXl^Xv = Id + 5i + QCv 

By Proposition EH we have that Cy e ^ i^^+^(f^i) for < A: < 

and hence Q/Iy : H^{n) H^+\n). Moreover, from (f02]) we have that 

QX^Xv = M + QF + Si + QR. 

We can then construct (Id + QF)~^, which is an elliptic pseudodifferential 
operator of order (its principal symbol is 1). Thus 

(Id + QF)-^QXpXv = Id + (Id + QFy^Si + (Id + QF)-^QR. 

We let 

(6.13) v = {ld + QF)-^QRf 

and note that v G H^^l). The other term (Id + QF)~^ Sif is microlocally 
smooth near (z,^). Thus, (z,0 ^ WF((Id + QFy^QX^Xyf) =^ (z,0 i 
WF(/ + v). Since 

WF((Id + QF)-^QXlrXvf) C WF((Id + QF)-^Q) n WF(Xf Xy/), 

the result follows. □ 

Remark 3. It is easy to see that the open dense sets of pairs (a, fc) that 
are dependent on I form a nested sequence. Moreover, one can eliminate 
the remainder function v by taking the intersection of all such open dense 
subsets. By the Baire category theorem, this limiting set of pairs will still 
be dense, but it is not clear if it remains open. 

Corollary 6.4. Suppose that Hi^T-j""*^ || < 1. Then in Theorem \2.^ we have 
v = 0. 

Proof. In this case the series ^"^oiKT^^y J converges to the identity plus a 
weakly singular integral operator F, which altogether is a pseudodifferential 
operator of order 0. In light of (j6.1ip this implies that Rq = 0, and hence 
R = 0. By (I6T2]1 and ([6T3]l we have v = 0. □ 
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Remark 4. Here we make a brief mention of how one might utihze the 
above results to detect a source / when partial data is known. Given z E fi, 
suppose that one can measure Xv4>z where 4>z = 6{x — z) is a point source 
centered at z. 

From the data Xyf, we may then compute XyXyfiz), since 
(6.14) X^Xvfiz) = {X^XvfAz)LHn,) = {Xvf,Xv<Pz)LHd^sn,,d^)- 

By Theorem 12.21 the graph {{z , Xy Xy f (z)) \ z E 17} provides an image 
which indicates, modulo a function v of some regularity depending on the 
regularity of a and k, the part of the wavefront set of / contained in A^'. 

However, in practice such a method of computing XyXy might be too 
computationally expensive. So in the case that ||rj~^iir|| is suitably small, 
one could try using a Neumann series truncated to one or two terms to 
approximate Xy. The author intends to consider this further in future 
work. 

Acknowledgements. This research was supported by NSF grant DMS- 
0838212. The results of this paper would not have been possible without 
some helpful and inspirational conversations with Gunther Uhlmann, and 
so for this the author expresses thanks. The author also thanks the journal 
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Appendix A. Relevant Results on Singular Integral Operators 



Suppose we have an integral operator of the form 
(A.l) 

Kf{x) = a{x)f{x) + / K{x, X - y)f{y) dy, K{x, x - y) = r~''4>{x, 



where / G ^/''(O), 6 = ^^Zy\ ; ™d r = \x — y\. The function (j){x, 0) is called 
the characteristic of the singular integral operator. We formally define the 
symbol ^>(x,^) of K by 

(A.2) ^(x,e) = 



e''''-^K{x,z)dz. 

It is easy to see by a change of variables that ^ is homogeneous of degree 
in ^. Letting oj = |||, we will write ^{x,ijj) from now on. It can be shown 
that if K[x, x — y) = (j){x, 9)r~'^, then 

1 



(A.3) 



^{x,uj) 



In 



cos 7 1 



ITT 

+ —sign (cos 7) 



d0 



where 7 is the angle between the vectors x and w. 

Consider the singular operator with a variable symbol, 
(A.4) 



{Af)ix) :=a{x)fix)+ 



^{x,e) 



^ \x — y 



nf(.y)dy 



e"-«<^A(x,a;)/(e)de, ^ 
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We introduce the class of those symbols that satisfy the condition 
(A.5) L>°$(x,w)G^r^(S"~^), Va:|a|</. 

Here the relation /3(x, a;)€-fr'(S"~^) means that 

(A.6) / \DZfi{x,bj)f <C, 0<|a|</. 

In this case, we say that /3(x,a;) belongs to H\§>"'~^) uniformly with respect 
to the parameter x. For symbols of singular integral operators that satisfy 
such a condition, we have the following useful theorem. 



Theorem A.l (Theorem XI.9.2, [13j). //«>a(x,w) G ^/,a where A > 
then the operator ( [^.^[ j is bounded in H^{W^). 



n~l 
2 ' 



For relating the characteristic (p{x,6) to its symbol ^{x,uj) we have the 
following theorem 

Theorem A. 2 (Theorem X.7.1, [13J). The symbol of a singular integral op- 
erator satisfies the relation ^{x,u)^H'^{S"'^^) if and only if the characteris- 
tic of this integral satisfies the condition cl){x, 0)€i?'(S"'~^) where I = X — ^. 

We also recall that the derivative of a weakly singular integral operator 
([S], IX §7) is given by 
(A.7) 



</>(x. 



^n— 1 



dy-f{x) / <l){x,e)ekdS{y). 



This formula holds for any / G L'^{Q) and for (p G (7^(17, S"~^). 

In dealing with weakly singular integral operators depending on a param- 
eter 9 and acting on functions / depending on x and 9\ it will be helpful to 
work with a particular type of space on which these operators work nicely. 
In particular, we will use expansions of functions in terms of spherical har- 
monics. Recall that any function g{x,9) G C°°(M" x S""-*^) can be expanded 
as a series 



(A.8) g{x,9) = Y,T.ci(^Hx)Yil{9), 

where 



m=0 k=l 



(2m -\- n — 2){m + n — Sy. 
[n — l)\m\ 

denotes the number of linearly independent spherical functions of order m. 
Furthermore, if g has compact support we claim that 

OO km,n 

m=0 k=l 
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In [22] it is stated that for g £ L°°(M"; C"+i(S"+i)) with compact support, 
we have that 

(A.IO) Yl llL-(R")II^Slli?i(S"-i) < °° 

m=0 k=l 

Under the assumption that g is compactly supported, we have that the L°° 
norm is comparable to the norm. Since all derivatives of g also satisfy 
(jA.lOp with llom^ ||loo(r") replaced by ||am^ Hia^ign-), we have that (|A.9P holds. 

Recah our definition of Uii^ x S""^) as the completion of C°°{n x S"^i) 
with respect to the norm || • ||-H;(r2xS"-i)i ^-iid C;(0 x §"~-^) as the completion 
of C°°{^ X S"^-*^) with respect to the norm || • ||c,(r2xS"-i)- 

The following proposition related to weakly singular integral operators 
and its applications in context will prove useful. 

Proposition A.l. Let A be the operator 

x-y 



with y, r, 9) compactly supported in x and y. Then for a constant C > 
depending only on n and I, 

a) //a e C72'+2(iR^ X M;; X X S^-i), then A : H\a) F'+1(M") is 
continuous with a norm not exceeding C\\a\\(j2i+2 . 

b) If a{x,y,r,9) = a' {x,y,r,9)(l){9) and also in C2'+2^ then 

II^IIh'(q)^h'+i(R") ^ c'II«'II(72*+2||<AIIhi(§"-i)- 

c) If f e niin) and a is as in (a), then A : Hii^ x S"-i) ^ H^+^i^) 
is continuous with ||^||-^;(f2xS"~i)— >-_H''+i(R") < C\\ 

d) Ifa = a(x, y, r, 9, rj) £ C°°(M5 x M^^ x x S^"^ x is compactly 
supported in x and y, then A : Tiiifl x S'"--^) V-i+iiW' x S"-i) is 
bounded with 

ll^llHi(CxS"-l)^-H; + i{R"xS"-l) < C'||a||c2;+2{K"X§"-1)- 

Proof. First note that in all cases, since a is compactly supported in x and 
y, we can let U d M" be such that a{x,y,r,9,r]) = for x,y ^ U. We 
can then replace a by aip without affecting the integral, where ip € C^(M) 
satisfies ip^r) = 1 for r € [0, diam(C/)]. Thus, we may as well assume that a 
is compactly supported in all variables. 

Consider the case that / is independent of 9 and a only depends on x and 
9. From the Calderon-Zygmund theory of singular operators, we know that 
for an integral operator K with singular kernel k{x, y) = (pix, 9)r~^ where 
we recall r = \x — y\, if (j) has mean value as a function of 9 for each x, 
then K is a well-defined operator on test functions, where the integration 
has to be understood in the principal value sense. Moreover, K extends to a 
bounded operator on satisfying ||i^|| < C sup^ ||</'(a;, •)llL2(§n-i) (Theorem 
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XI. 3.1, [H]). As a remark, the extension can be considered as a convolution 
in the sense of distributions, and then need not have mean value in 0. 

Let (ji, j2) • • • 1 Ji+i) be a multi index. To make notation a bit more con- 
sistent, let ao = ot. Consider the derivative dx^^A, which by ()A.7j) and 
([22], Proposition 1) consists of a bounded term ai{x)f{x) plus the integral 
operator with kernel 

/A 11^ a a(x, g) _ (1 - n)gj,a + de^^ a ^ d^^^ a{x, 9) 

(A.lij (Jx. , - — h 



Jl j,n—l j,n j,n—l 

Letting (j)i{x,6) := (1 — n)9j^a{x,9) + dg._^a{x,9), since a G C'2«+2 ^^^^ -g 
compactly supported in x, we have that the symbol 



<^i{x,u) = J e--"" j^^J dz 

of ^1 belongs to C^'^^(M" x S"^-*^). Since i?i>i is compactly supported in 
X, a2(/'iG^^'+^"l^l(S""^) C -H'^+^-t(S"-i) for ah < I7I < / and for 
some fixed A > (in particular we could take A = ^). By Theorem 
IA.21 $i(x,a;) G By Theorem lA.ll we have that the integral kernel 

(j)i{x , 9)r~'^ corresponds to a singular integral operator that is bounded on 
HK For the second term in (lA.lip . which is a weakly singular integral ker- 
nel, we have that ai{x,9) := dxj_^a{x,9) € C^'^^. Similarly as before, we 

compute dx^^ ^ ^ which corresponds to an operator with a bounded 

multiplier 02 (x), a singular integral operator, and a weakly singular integral 
operator. It can be shown analogously that the symbol <I>2(2;,^) correspond- 
ing to the characteristic 02 (x,^) of the singular integral term belongs to 
£?i-i^x. Thus modulo a weakly singular integral operator, the operator A2 
with kernel (j)2{x,9)r~^ is bounded on H^~^. We then focus our attention 
on the weakly singular integral operator that remains. 

After repeating this process a total of / -|- 1 times, which involves / -|- 1 
differentiations, the remaining weakly singular integral operator has a kernel 
a;+i(x, 9)r~^^^ with a/+i(x, 9) G C^^^. We can then proceed as in the proof 
of ([22j. Proposition 1) to obtain that this term is bounded on L^(J7). In 
particular, we use the criterion from Calderon Zygmund Theory which states 
that if K is an integral operator with integral kernel A;(x, y) satisfying 

(A.12) sup \k{x,y)\dy<M, sup \k{x,y)\ dx < M, 

X J y J 

then K is bounded in with a norm not exceeding M (|23j. Prop. A. 5.1). 

Now we want to bound the operator norm ||^||H'(r2)^_f/'+i(R") this sim- 
pler case. Let <j)i be the characteristic of the ith. singular integral operator 
obtained by the above process with symbol <I>j and Ai the corresponding 
singular integral operator. Note that Ai is bounded on H^~^^^(ft). One can 
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explicitly compute that 

(A. 13) Oi = dx^^xj^-xj^a, 

(A.14) (pi = {l- n)ej^ai-i + de^ai^i, 

(A.15) a,{x) = j Qi_i(x, 6)6 dS{y), for 1 < i < / + 1. 

Also define 



(A.16) R^f{x)■.= I ^^^f{y)dy. 



We have 

||<Ai|lc"+i-' < C'll"llc"+i 
||ai||ci+i-i < C||a|lc"+i 
(A.17) ||ai||ci+i-» < C\\ai-l\\cl+l-^ < C\\ai-i\\ci+i-(i-i) < C'\\a\\ci+i 

If = 1> then ([I3j, Thm XI. 3.2, Thm XI.9.2) imply that for 1 < i < 



Mi/ll_f/'-»+i{R") - X] ll-^f^i/llL2(IR") 
|/3|</-i+l 

<^ E E^^Pll^-''''^*(^'-)ll^^(S"-)ll/ll/^Kn) 

|/3|<i-i+l7</3 

< C||(?!>j||f;i-^i+l 

(A.18) <C\\a\\ci+i. 



AN INVERSE SOURCE PROBLEM WITH PARTIAL DATA 23 

Now we estimate using ()A.18h with the understanding that we 

sum over all indices ji, j2, • • •• Again, assume that = 1- 

\\Af\\Hi+. = \\AfU2 + \\d,^^Af\\Hi 

< \\a\\co + ||ai/||^^! + pi/ll/fi + WRifWrn 

< \\a\\co + llaillc-i + ||a||c'+i + II^i/IIl2 + \\dx,^Rif\\w-i 

< ||a||c'+i + ll«llco + ll"i|lco + llaillc + ll«2/lbi-i + \\A2f\\Hi-i + 11^2/1^;- 

< ||a||ci+i + l|a||co + ll"i||co + Ikillc' + ll«2||c'-i + ll^2/||iyi-i 



i+i 



^ ll-Ri+l||L2 + !l"llc'+i + ^ Il"i-l|lc0 + ll«i|lc'+i-' 



i=l 
l+l 



^ ll"«+lllco + ll"llci+i + ll«llc'+i 



< \\a 



i=l 



Thus in the simplified case where a only depends on x and 9 and / is 
independent of 0, we have 

(A.19) ll-4||//!(f7)^//!+i(Mn) < C'||a||c'+i- 

To extend to a = ag depending also on y and r, we use a first order 
Taylor expansion in y and r centered at y = x and r = 0, similarly to in 
([22], Prop. 1), to get 



ao{x,y,r,9) = ao{x,x,0,9) + ^ r^'^^iy - x^ d^d^ao{x, x + t{y - x),tr,9) dt 

l/3| + l7l=l ° 

= ao{x,x,0,9)+ rl'^l(-re)T / di!d^ao{x,x + t{y - x),tr,e) dt 

l/3| + l7l=l ° 

= ao{x,x,0,9)+r V (-1)1^1^^ [ d^d^aoix,x + tiy - x),tr,9) dt. 



-|7l=l 

We can then write 

(A.20) ao(x, y, r, 9) = ao{x, x, 0, 0) + r-fi{x, y, r, 9) 

where 71 G C^'"^^. After dividing by r"~^, the first term in ()A.20p maps 
to H^~^^ by the previous argument. The second term corresponds to 
an integral operator with kernel 7i(x, y, r, 0)r~"+-^. If we differentiate this 
with respect to x, we get a weakly singular integral operator with kernel 
ai(x, y, r, 0)r^"^"'^ where qi € C^'. Now repeat as before, writing 

(A.21) ai(x, y, r, 9) = qi(x, x, 0, 6) + r72(x, y, r, 9) 
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where 72 G C'^''^^. The first term ai{x,x,0,0) corresponds to a bounded 
operator Ai : H'-^^{i}) — )• H\W^). Moreover, j2{x,y,r,9)r~^~^'^ can be 
differentiated with respect to x to obtain a weakly singular integral operator 
with kernel a2{x,y,r,9)r~"'^^ where 02 £ C^'~^. 

After repeating this process a total of I times, we have a remainder term 
that is a weakly singular integral operator with kernel ai{x,y,r,9)r^^^^ 
where a; G C^. Write 

ai{x, y, r, 6) = ai{x, x, 0, 9) + rji+iix, y, r, (9) 

with 7i+i G C^. Then 7^+1 corresponds to the operator 7;_|_i(x, y, r, 6')r~""'"^, 
which we can differentiate with respect to x to obtain a weakly singular 
operator that is bounded on with a bound not exceeding ||7/_|_i||f7i by 
using the estimates in (|A.12p and applying the Calderon Zygmund theorem. 
Since each weakly singular integral operator with kernel aj{x, x,0,9)r~''^~^^ 
is a bounded map from H^~^{Vl) — )■ i?'~-'"^^(M"), we combine the remainder 
terms together to get that A : H\n) i7'+^(M"). 

More explicitly, let Ai be the weakly singular integral operator with kernel 
ai{x, X, 0, 9)r~''^~^^ and Ri the integral operator with kernel 7j+i(x, y, r, 6')r~""'"^. 
In particular, Ai = dx Ri-i- We will also need the straightforward estimates 

(A. 22) ||aj||c-i-i+i < ||a||ci+i+i, ||7j||c"" ^ i|a||cm+2i-i . 

For 11/11//! = 1 we have 

P/lbm<Po/||Hm + ||i?o/lbm 

< ||q!oIIc'+i + hilico + WMfWm + WRifWm 

< ||ao|lc"+i + hilico + llaillc + II-Ri/IIl2 + Wdx^^RifWrn-^ 

< ||ao||c'+i + ll«i|lc' + IItiIIc + II72IIC0 + WMfWm-^ + 11-^2/11^;-! 



ai\\ci+i-i + ||7j+i||(70 

i=0 

l-l 

< \\ai\\ci + ||^i/||L2 + \\dx,^_^^Rif\\L^ + \\ai\\ci+i-^ + ||7i+ilbo 

1=0 

l-l 

< ||a||c.2i+i + ||7z+i||co + ||7i+i||ci + ^ + ||a||(^2i+i 

i=0 

l-l 

~ ||o||(72i+i + ||a||(72i+2 + ^ ^ I|q^IIc'+i+* ~^ llo^llc^'+i 

i=0 

< ||a||c2i+2. 
This proves (a). 
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Now consider if a(x,y,r,9) = a'{x,y, r,6)(l){6). Then 

(A.23) (1 - n)eja + d0.a = {l- njOja'^) + a'de.4> + (fide. a' . 

In short, for each term in the decomposition of A by differentiation, (j) is 
differentiated exactly once. Therefore 

(A.24) ll^l|j?((f2)-).i7'+l(]Rn) < C||a'||c.2i+2||<;/!)|ii-l(Sn-l), 

which proves (b). 

For f eniinx S"-^) depending on 9 as well, we expand / as a series 

m=0 k=l 



Then 



[Am 



a{x,y,r,e) 

— 1 f{y,d)dy 



m=0 k=l 



EE/ 

m=0 ifc=l 

Hence 

m/iiif'+i(R"x§"-i) < E E 



a{x, y, r,e) 



r 

CO km,n /. / 



r-n— 1 



m=0 

OO km.n 



a{x,y,r, e)Y^l{0) (k) 



r.n—1 



aZKy)dy 



i?'+i(M"xS"-i) 



m=0 it=l 

OO ^m.n 

< C\\a\\c2i+2 ^ ^ \\^m)i\\Hi{S»--L)\\a^m\\Hi'{^) 

m=0 k=l 

< C\\a\\c2l+2\\f\\■u^(^a)■ 

Here the operator A^^n is given by [A^lgjix) = J ^^^^^^^^0^^ g{y) dy . 
This proves (c) . 

Finally, if a = a{x, y, r, 0, rj) is C°° with compact support, we can expand 
it as a series 

OO km.n 

(A.25) a(x,y,r,0,r?) = ^ ^ 6«(x,2/,r,0)y^^) (77). 

m=0 k=l 
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Note that a G Cj{W' x R" 
/ G UiiVL X S'^-i), we have 

ll^/l|-Hi+i(IR"xS"-i) 

= X] ll^ifi^,nllHi(§"-i) 

mi=0 fci=l 

oo ^m-i,n 

= X] ll^ifi^,nllHi(§"-l) 

mi=0 fci=l 



X M X S"-^ X S"-^) for ah j > 0. Then for 



bl^l\x,y,r,e) 

z^i fiy,0)dy 



_ff'+i(IR"x§"-i) 



oo km2,n Akl) 



m2=0 k2=l 



bt,\x,y,r,d)Yi'2i{e) 



y.n—1 



_fi''+i(R"x§"-i) 



OO ,ri OO ,^1- 

mi=0 fci=l r?i2=0 A;2=l 



□ 



«llC2,+2ll/IIW,(f^xS"-l)- 
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